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\ Abstract. The location of the images in a multiple-image gravitational 

\^ . lens system are strongly dependent on the orientation angle of the mass 

distribution. As such, we can use the location of the images and the 
photometric properties of the visible matter to constrain the properties 
of the dark halo. We apply this to the optical Einstein Ring system 0047- 
2808 and find that the dark halo is almost spherical and is aligned in the 
Q"^ ' same direction as the stars to within a few degrees. 
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Q . 1. Introduction 

^ , 

Q^, Numerical simulations of Cold Dark Matter (CDM) have been very successful in 

reproducing the observed large scale structure of the universe. The CDM model 
predicts that the dark matter (DM) haloes of today's galaxies are assembled 
through successive mergers of smaller haloes. Simulations using only dark matter 
53 . predict that the haloes should be quite prolate, however it is not clear how gas 

and/or stars interacting with the dark matter will change the shape of the halo. 
Studies have suggested that the DM halo can become more or less cuspy (El- 
^ I Zant, Shlosman, & Hoffman, 2001; Tissera & Dominguez-Tenreiro, 1998) and 

. rounder (Evrard, Summers, & Davis, 1994; Dubinski, 1994) after the interaction 

with stars and gas. An important test of galaxy formation and evolution models 
will be to compare the shape and profile of galaxy haloes with observed haloes. 
Thus, simple questions such as: "Do we expect the visible and dark matter to 
be aligned in elliptical galaxies?" and "Is the dark matter density in the central 
regions changed by the gravitational dominance of the stars?" must be answered 
with observations. For instance: the Milky Way, despite being a spiral galaxy, 
appears to have an almost spherical halo (Ibata et al., 2001). 

Gravitational lensing offers a method to tightly constrain the shape of DM 
haloes in the population of medium redshift (0.1 < z < 1.0) lens galaxies. The 
image positions in a lens system are highly sensitive to the orientation of the 
overall mass profile. Keeton, Kochanek, & Falco (1998) showed that the overall 
mass distribution is typically aligned with the visible matter using a sample of 
lens galaxies and a simple SIE mass model. However, depending on the lens 
galaxy, the stellar mass can contribute a substantial fraction of the total mass 
inside the image. The extreme case is the lensed QSO 2237-1-0305 where the 
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dark matter constitutes only 4% of the projected mass inside the images (Trott 
& Webster, 2002). In this case we expect the visible matter orientation and the 
total matter orientation derived from a lensing analysis to be very similar. The 

logical next step is to use a more complicated (stars + halo) model for the lens 
galaxy to determine the properties of the DM halo alone. 

In this paper we use an implementation of the LensMEM algorithm (Walling- 
ton, Kochanek, & Narayan, 1996) and a stars+halo lens model to study the 
optical Einstein Ring 0047-2808 (Warren et al., 1999, 1996) using data from 
the HST. This system is well suited for the study because it is an isolated lens 
galaxy so we expect any external shear contributions to be small. The system 
is a z = 0.485 elliptical which is lensing a background starbursting galaxy at 
z = 3.6. 

The algorithm we employ performs a non-parametric source reconstruction 

to match the observed data for a given lens model. The goodness-of-fit of the 
model is calculated using a taking into account the degrees of freedom used 
in the source. In this paper we assume B.q = 70 kms~^Mpc~^ and {^m, ^x) = 
(0.3,0.7). 



2. Method 

The data was reduced as described in Wayth et al. (2002). The final image of 
the "ring" is 133x133 with 0.05" pixels as shown in Figure 1. The lens galaxy 
was best fit with a Sersic profile, where the surface brightness as a function of 
radius r is E = S1/2 exp{— B(n)[(r/ri/2)^^'* — !]}• The parameter n quantifies 
the shape of the profile: the values n = 0.5, n = 1, and n = 4 correspond to 
the Gaussian, exponential, and de Vaucouleurs profiles. Profiles with larger n 
are more cuspy. B{n) is a constant for a particular n and we used the series 
asymptotic solution for B{n) provided by Ciotti & Berlin (1999). Additional 
parameters used for the light profile are the axis ratio (q) and orientation angle 
{6s)- The fitted parameters are shown in Table 1. 

Ri/2 (pixels) 21.69 

Ei/2 (counts) 0.7 

q 0.693 

0s n 125 

n 3.115 

Table 1. Photometric parameters for the lens galaxy. 

We model the galaxy stellar component with fixed parameters from the 
photometry and allow only the M/L to vary. The halo is modelled as a Pseudo- 
Isothermal Elliptic Potential (PIEP) with a finite core. The PIEP model is 
defined by the lensing potential V' = b[r^ + (1 — e)a;^ -|- (1 -|- e)y^]^/^ where rc 
is the core radius, b is the mass scale (Einstein radius) and e is the cllipticity. 
An additional parameter is used for the orientation angle {9h, measured anti- 
clockwise from horizontal). It is worth noting that the lens can be fit with the 
PIEP model alone (without a core) with the parameters b = 1.165, e = 0.08 and 
Oh = 129. We use this mass scale for the halo model. The source in this system 
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Figure 1. Data and model for the 0047-2808 lens. Left: the HST 
image and image with the galaxy subtracted. Right: Model image and 
reconstructed source for the PIEP model. 

actually has two distinct components. The two-component model explains the 
location and brightness of all features in the image with a standard lens model. 
Figure 1 shows the model source and corresponding image for the plain PIEP 
model. 

The mass enclosed inside the image is tightly constrained by the Einstein 
radius. We use this constraint to normalise the stellar M/L for a halo of a given 
core radius. A large core is equivalent to a constant M/L mass model, whereas 
a small core will generate an unrealistically low M/L for the observed stellar 
component of the lens. 

In preliminary tests, we found that we cannot fit the data for rc ^ 7" (42kpc 
physical scale length) i.e. constant M/L models cannot fit the data. Therefore 
we have restricted our analysis to halo core radii < 7". For the range of allowed 
core radius values, we have calculated the range of halo ellipticity and orientation 
angle which produce acceptable fits to the data. 

3. Results 

Figure 2 plots the acceptable range of halo ellipticity and orientation angle as a 
function of core radius. On the left, we see that the halo ellipticity is consistently 
less than the stellar ellipticity. For 1.5" < Tc < 2.5", the data permit a halo with 
projected mass density which is circular, although in all cases the best solution 
has a halo with non-zero ellipticity. 

On the right of Figure 2, the plot shows that the halo orientation angle 
is independent of the core radius and is in the same direction as the projected 
stellar major axis (within errors). The acceptable range of orientation angles for 
1.5" < Tc < 2.5" are for non-zero ellipticity. 

4. Conclusion 

By using a lens model which separates the stars from the halo, we have been 
able to determine some of the basic properties of the dark matter halo in the 
lens system 0047-2808. We find that the projected ellipticity of the halo is not 
circular, but is substantially rounder than the observed stellar ellipticity. A 
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FigLirc 2. Constraints on the dark matter halo. Asterisks indicate 
the best fitting parameter values. The error bars indicate the range of 
the parameter which will produce acceptable (la) values. Left: The 
acceptable range of ellipticity (e/^) for various core radii. Right: the 
acceptable range of orientation angle {Qh) 

small range of halo core radii values (1.5" < rc < 2.5") allow the projected halo 
mass to be circular. 

The halo's core, modelled as a constant density region, must be < 7" to fit 
the observation. The core size could be further constrained by applying realistic 
limits to the stellar M/L which we intend to do in further work. 

Finally, we find that although the halo is less elliptical than the stars, the 
orientation angle of the star's and halo's major axis are the same within errors. 
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